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Abstract 

Radiative corrections to the muon decay spectrum due to soft and virtual electron- 
positron pairs are calculated. 
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1 Introduction 

The experiment TWXST [1,2] is currently running at Canada's National Laboratory TRI- 
UMF. It is going to measure the muon decay spectrum [3,4] with the accuracy level of about 
1 • 10~ 4 . That will make a serious test of the space-time structure of the weak interaction. 
The experiment is able to put stringent limits on a bunch of parameters in models beyond 
the Standard Model (SM), e.g., on the mass and the mixing angle of a possible right- 
handed VT-boson. To confront the experimental results with the SM, adequately accurate 
theoretical predictions should be provided. This requires to calculate radiative corrections 
within the perturbative Quantum Electrodynamics (QED). Here we will present analytical 
results for two specific contributions related to radiation of virtual and soft real electron- 
positron pairs. The corrections under consideration are of the order O (a 2 ), where a is the 
fine structure constant. 

The contributions of virtual u, + ur , r + r~, and hadronic pairs were found [5] to be small 
compared with the 1 • 10~ 4 precision tag of the modern experiments. The contribution of 
e + e~ pairs is enhanced by powers of the large logarithm L = ln(m^/mg) ~ 10.66. Analysis 
of the leading and next-to-leading terms from this correction in Refs. [6,7] has shown that 
the numerical effect is not as small as for other leptonic flavors, and it should be taken 
into account. Comparison of the leading and next-to-leading contributions revealed a poor 

1 A certain part of this work was performed in University of Alberta, Edmonton, Canada 



Preprint submitted to Elsevier Science 



1 February 2008 



convergence of the series in L. Calculation of the terms without the large logarithm was 
found to be desirable. 



Within the Standard Model, the differential distribution of electrons (summed over electron 
spin states) in the polarized muon decay can be represented as 



= T [F(x)±cP,G(x)], T = 



- 



dxdc ^ V-V — M-V~VJ> ^ 1927r3 ' 

c = cos6>, x = —5- -, x <x<l, x = ^- J f — -, (1) 

+ m z e + m z e 

where m M and m e are the muon and electron masses; Gf is the Fermi coupling constant; 9 is 

— # 

the angle between the muon polarization vector P M and the electron (or positron) momen- 
tum; E e and x are the energy and the energy fraction of e ± . Here we adopt the definition 
of the Fermi coupling constant following Ref. [8]. Functions F(x) and G(x) describe the 
isotropic and anisotropic parts of the spectrum, respectively. Within perturbative QED, 
they can be expanded in series in a: 



F(X) = /BornOr) + ^h(x) + (^M*) + + O (<**) , (2) 

and in the same way for G(x). Among different contributions into the functions F(x) and 
G(x) (see Ref. [6] for details), there are ones related to the production of electron-positron 
pairs. In this Letter we will consider the effect of soft and virtual e + e~ pairs. 



2 Soft e+e Pairs 

The process of real pair production doesn't reveal any infrared singularity, contrary to the 
case of photon radiation. Nevertheless, a separate consideration of soft pair emission can be 
of interest. In fact, e + e~ pairs with energy below a certain threshold can't be observed in 
experiments with muons decaying at rest. So, the corresponding contribution is a specific 
correction to the measured decay spectrum. Moreover, the behavior of the real pair emission 
in the soft limit is not smooth. An integration over the domain between the threshold of 
real pair production and a certain cut on the maximal energy of the soft pair is desirable. 

The maximal energy of the soft pair is defined by the parameter A, which is assumed to 
be large compared with the electron mass: 

^pair < « A « 1. 3 

2 m M 

Due to the smallness of the pair component energies, the matrix element M of the process 
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H (p) — >e (q) + v^n) + u e (r 2 ) + e + (p+) + e (p_) (4) 

can be expressed as a product of the matrix element M of the hard sub-process (the 
non-radiative muon decay) and the classic accompanying radiation factor: 



M = M —v( P+ )Yu(p-)J^ k = p + +p-, (5) 

where p+- are the momenta of the positron and electron from the created pair. The radi- 
ation factor reads 



j = 2a la (q) 



Performing the covariant integration of the summed over spin states modulus of the matrix 
element over the pair components momenta, we obtain 



E \v{ P+ )l»u{p_)\ 2 = A(p^_+p^_ - ^ gT), 
spin 

J P+P- 4 



It is convenient to parameterize the phase volume of the total pair momentum as 

d 4 k = dk k 2 d\k\dQ k = 7rdk dk 2 \JkQ — k 2 dc k d^ k , (8) 

where a trivial integration over the azimuthal angle can be performed: / dip k — * 2tt. Now I 
integrate over the total pair momentum with the condition (3) (k = E pair ). In this way I 
got the following result for the soft pair contribution: 



dr SP jpBorn dr Born 



—. — j — = To [fo(x) ± cP^go(x)} + O 



m 



2 ■ 



dc dx dc dx dc dx ' ' \ ///;, 

fo(x) = x 2 (3-2x), g (x)=x 2 (l-2x), 
5 s P a 2 r 1 , 3 . 2, . , , /61 _A 223 



— \n 3 A--\n 2 A + \nA(— - ((2)) - 
12 3 U8 V V 27 



3tt 2 

+ |C(2) + 2C(3) J 

oo i 2 

lnA = L + 21nA, CM = E C(2) = y 



(9) 
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So we calculated explicitly all the terms in <5 SP except the ones suppressed by the small 
factors (a/7r) 2 m 2 /m 2 and (a/7r) 2 A. 



3 Virtual e + e Pair 



We will use here the substitution suggested by J. Schwinger for the photon propagator 
(with 4-momentum k) corrected by a one-loop vacuum polarization insertion: 



I .a}dv<f>(v) 1 M 2 = ^±, (10) 



k 2 -X 2 +iO 7T J 1-v 2 k 2 -M 2 + iO' l-v< 


Hv) = l-l(l-v 2 )(2-v 2 ), 
where m 2 is the mass of the fermion in the loop. 

The standard technique of integration over Feynman parameters can be used here. We are 
interested in the region of electron energy fractions z ^> m e /m^. Analytical expressions 
for the relevant integrals in this region are given in Appendix A. As concerning the region 
of small electron energy fractions (z ~ m e /m M ), it requires a more accurate treatment. 
But the differential width there is rapidly decreasing (see i.e. the Born-level functions in 
Eq. (9)), and the contribution of this region into the total width is suppressed by the mass 
ratio squared. 

Formally, we have an ultraviolet singularity in the virtual pair correction. The Fermi theory 
is not renormalizable in the general case. But for the muon decay everything is safe, since 
the standard renormalization of the electron and muon wave functions removes the singu- 
larity [9]. Note that we need to use here only the pair contribution into the renormalization 
constants. They can be found easily from the calculation of the virtual pair corrections to 
the ee and \x\i vertexes (see Appendix B), where we had 



A 2 A 2 
In — — -2\4(0), In > -2V,(0), 



V,(0) = -U? + f 6 L - H - I C (2), V.(0) = + | C (2), (ID 

where A is the ultraviolet cut-off; see Appendix A for quantities D ejfl . For the muon decay 
we use half a sum of the above substitutions: 



5 ln ! + 5 ln f^-^ 0) - v;(0) - (12) 

The same logic was applied for renormalization of one-loop corrections to muon decay [14]. 
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I got the following result for the virtual e + e pair contribution: 



pVP 

dc dx 
where 



= r r 



ot 



2tt) 



ti%\x)±cP»g^ , {x) + 



(e+e"). 



(13) 



/iwrt V) = /o(a:)W(x) - 2x 2 lnxL - 2x 2 In" x - 2x 2 Li 2 (1 - x) 



2l„2 



+7x lnxH — xlnxH — In a; — 



3 ' 3 3(1 -x) 

2 o, , 2 



lnx, 



#2 



virt 



^x) = g (x)W(x) — -x 2 lnxL — -x 2 ( Li 2 (1 — x) + In 2 x 
+ — ^x 2 lnx — -xlnx — - lnx + 
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2 
3' 



2 
3 



lnx, 



3(1 -x) 



X 



oo 4 3Q7\ 8 4 38 

+ - lnx - -C(2) - — JL - -S 1>2 (1 - x) + -Li 3 (1 - x) + y Li 2 (1 - x) 



38 



265 



-In x lnxLi 2 (l — x)H In x £(2) lnx —lnx 



9 3 
4 22 517 

+ ^C(3) + #C(2) + " 1, 
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27 



Li 2 (x) = - / dy 



ln(l-y) 
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Li 3 (x) = dy 



27 



(14) 



L12 (y) 



q / x _ 1 / , In (1 -y) 

o y 

It is worth to note that the sub-leading virtual corrections don't factorize before the Born 
functions /o(x) and go(x). 

By integration over the energy fraction and the angle we receive the corresponding contri- 
bution to the total muon width: 



.VP 



dx 



w 

-1 o 

20063 61 



pVP 



dc dx 



r 



16 



a 
2~ix 



.9 3 V 36 3 SV V 
-5.0497- 10~ 5 T . 



(15) 



This quantity was calculated earlier in Ref. [10] by numerical integration using dispersion 
relations: 



r vp ([10]) « -5.1326- 10" 5 r , 



(16) 
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which is close but different from my number. The reason for this discrepancy will be in- 
vestigated elsewhere. At least part of the difference can be due to terms proportional to 
(a/7i) 2 (ml/mf l )L n , which were omitted in my calculation. 

The correction to the forward-backward asymmetry of the decay can be found also: 



F VP 
1 FB 



i o 
dc — J dc 
n -l 



dx 



pVP 



+ i C (2))L-^4c(2)-SC(3) 



7776 



dc dx 

-C(2) 
2 SV ; 9 



1 



13 



/647 



54 54 ' V648 
-1.17- 10" 5 r . 



(17) 



4 Numerical Results and Conclusions 



The relative effect of the soft pair correction depends only on the cut value. It is shown in 
Fig. 1. The soft pair approximation (3) is not valid for values of A close to the threshold of 




Fig. 1. The relative effect of soft pair corrections versus the cut value. 

real pair production and for large A ~ 1. But it can be used there as a simple estimate. So, 
by taking A = 1 we make an estimate of the order of magnitude of the total contribution 
due to real e + e~ pairs (here the estimate is about two times the true value). For very small 
values of A the correction should vanish in any case, so the approximation is really safe 
there. 

Let us define the relative contribution of the virtual e + e~ pair corrections in the form 



xVTV \ ( a \ 2 /2,virt ^ ( X ) + C-f/^lvirt ' ( x ) , 

6 W = fcJ /„M + eP„ JoW ' (18) 

The dependence of this function on the electron energy fraction is shown in Fig. 2 in different 
approximations for P M = 1, c = 1. The dependence on c is very weak, because the main 
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part of the correction is factorized before the Born-level functions and cancels out in the 
ratio. The leading logarithmic (LL) approximation takes into account only the terms of the 
order O (a 2 L 3 , a 2 L 2 ) , the next-to-leading logarithmic (NLL) approximation includes also 
the O (a 2 /. 1 ) terms, and the next-to-next-to-leading approximation (NNL) represents the 
complete result. 
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Fig. 2. The relative effect of virtual pair corrections versus electron energy fraction in different 
approximations . 



The third power of the large logarithm cancels out in the sum of the virtual and soft pair 
contributions: 



pSVP 

dc dx 



r °(0 



,(e+e- 
i2,SV 



(X 



x) + 



m 



m 



A 



(19) 



where 



fvsv \x) = fo(x)U(x) - 2x 2 InxL - 2x 2 In 2 x - 2x 2 Li 2 (1 - x) 

2 2 2 2 

-InrrH — x\nx + 7x lnxH — lnx, 

3(1 — x) 3 3 

fl&v ) (x) = go (x) U(x) - - x 2 In xL - -x 2 In 2 x - -x 2 Li 2 ( 1 - x) 

2 , 2 , 13 2 , 2 , 

H — ; -mx xmx-\ x mx mx, 

3(1 -x) 3 9 3 

/l 2 2 \ /4 32 

U(x) = f - + - In A - - \nx\L 2 + ( - In 2 A - — In A 

- -Li 2 (1 - x) - - In 2 x + — In x - — - -C(2) )l 

3 V ; 3 9 6 3 SV V 

+ ^ In 3 A - - In 2 A - -C(2) In A + — In A 
9 9 3 SV ; 27 

4 8 8 8 

+ -Li 3 (l-x) - -Si )2 (l -x) - -ln 3 x- -lnxLi 2 (l - x) 
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38 T . , . 38 . o 
+ —U 2 (l-x) + — \n 2 x 

+ ^ + 6C(2) + 4C(3). 




(20) 



I checked that the leading and next-to-leading terms in the above formula agree with 
the corresponding contribution obtained within the fragmentation function formalism in 



In this way we simulate the experimental set-up with a certain energy threshold for regis- 
tration of pairs, while events with pair production above the threshold (with several visible 
charged particles in the final state) are rejected. 

If the radiation of real pairs is completely forbidden by kinematics (or experimental con- 
ditions), only the virtual corrections (14) contribute. That happens, for instance at large 
values of x > 0.99. 

Thus, two contributions to the total set of radiative corrections for the muon decay spectrum 
are presented. They are required to reach the level of the theoretical accuracy below 1 • 10 -4 . 
The formulae can be used for semi-analytical estimates and as a part of a Monte Carlo 
code to describe the pair production contribution to the decay spectrum. The formulae are 
valid also for pair corrections to leptonic r-decays. 
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Appendix A 

List of integrals for the virtual pair correction 

Here we give the list of integrals over Feynman parameters, which are relevant for the 
calculation of the virtual pair correction to the muon decay spectrum (z ^> m e /m^ is 
assumed) . 



Refs. [6,7]. 



<-> 

D 



= I [y 2 l I - T2 L * + ( Li2 (1 " z) + T + c(2) ) Lz " Lis (1 " z) 

+ 2Si, 2 (1 - z) - hi 2 {l-z) + \- y C(2) - C(3), " , 
= l - \^L\ - h z + 2Li 2 (1 - z) + 2 + 4C(2)] , 




( 



yxB 




— In zL z — Li 2 (1 — z) + In 2 z + - In z , 



D 



— z 



8 



V D ' 

, y 2 xB 

[ D ' 

, y 2 x 2 B 

{ D ' 



A* 



17 13 

i ^- T ^ + Li 2 (l-z)- y + 5C(2) 



3(1 -z) 

z 

3(1^ 



1 1 17' 

--lri2;L z - -Li 2 (1 - z) + - In 2 z + — \nz 

^-^ + ^zlnz^L z + ^U 2 {l-z)-^ln 2 z 



-^lnz-lnz- y(l-z) 



1 


1-2 


3(1 -z) 


[ 4 



19 



+ ^- \nz - \nz - -(1 - z)jL, - -Li 2 (1 - z 

1 Q 10 

| In 2 ; + In 2 ^ + j|(2 - z) hi z - -£(1 - z) + 5(1 - 2 )C(2) 

-^[(,-1- £1^-^(1-,) + fin 2 , 

5 , ol 19,., vl 

zlnz + 2\nz H (1 - z) 



12 

L + 2 In z. 
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(A.l) 



I used above a short notation for the integral over three Feynman variables: 



dv (j){v) 



(F(y,x,y)) = J ^— J J ydy J dxF(v,x,y), 



o o 



and 



(A.2) 



D = y 2 P 2 + (l-y)M 2 , 



x 



'■ml + (1 - x) 2 m 2 e + Bx(l - x), 



D, = y 2 ml + (1 - y)M 2 , D e = yW e + (1 - y)M 2 , 



M 2 = ^, <P(v) = 2 -- 1 -(l-v 2 )(2-v 2 ). 



(A.3) 



Appendix B 

Asymptotic expressions for the muon form factor 



Using the Schwinger substitution (10), I reproduced the known [12,13] asymptotic expres- 
sions for the O (a 2 ) virtual pair contributions into the Dirac form factor of muon: 



Fi 4 > a \ mi ,m 2 ,Q 2 
'265 



216 



mi,m 2 <CQ 2 



7T 



36 



19 
72 J 



-L 



(B.l) 
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£>(!) = — --C(2), 

, 3355 19, , 1,, x 
= 1296 + 36 C(2 > - 3« 3) ' 




67 _7_ 

54 ~ 36 




C(3), 



where m\ = is the muon mass; rri2 is the mass of the fermion in the loop; e and e 2 is the 
muon and fermion charges, respectively; — Q 2 is the square of the momentum transferred 
in the spacelike region: — Q 2 = (pi — p 2 ) 2 < 0, where pi and p 2 are the initial and the final 
muon four-momenta. 
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